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ABSTRACT 

We present results from a study of ab-type RR Lyrae variables in the Large Magellanic 
Cloud using the extensive dataset from phase III of the Optical Gravitational Lensing 
Experiment (OGLE). The metallicities of the RR Lyraes, determined from the periods 
and amplitudes of their light curves, reveal a statistically significant radial abundance 
gradient that is approximately one-half of what is seen in the disks of the Milky Way 
and M33. The RR Lyrae abundance gradient agrees with that of the old and metal- 
poor LMC globular clusters. The reddenings of the OGLE RR Lyraes have been 
calculated using their minimum light colors and reveal a mean value of E(V-I) — 0.12 
± 0.05, where the quoted uncertainty represents the standard deviation of the mean. 
The distribution of RR Lyrae extinctions across the face of the LMC is well-correlated 
with the distribution and emission intensity of CO clouds based on recent millimeter 
wave surveys. In addition, we find that the old LMC globulars tend to be located in 
regions of low extinction. This underscores the need to survey the higher extinction 
regions with the specific aim of increasing the sample of old LMC globular clusters. 
Finally, we examine the distance distribution of the RR Lyraes in order to probe the 
structure of the LMC and investigate the possibility that some of the RR Lyraes may 
reside in a kinematically hot halo population. In addition to calculating a mean LMC 
distance of (m — M)o = 18.55 ± 0.10 mag, we conclude that some fraction of the RR 
Lyraes in our sample could be members of the LMC halo. 

Key words: stellar populations - RR Lyrae - LMC - dwarf galaxies. 



1 INTRODUCTION 

The ancient stellar population of the Large Magellanic 
Cloud (LMC) is one of the keys to unlocking the formation 
and early evolution of this very important dwarf irregular 
galaxy. For decades, the metal-poor globular clusters in the 
LMC have been used to study its early history because they 
represent distinct points in space and time (van den Bergh 
2004; Colucci et al. 2011). Knowledge of the ages, chemi- 
cal abundances, and kinematical properties of these clusters 
has led to a great many advances in our understanding of 
the LMC's early evolution. However, based on the latest cen- 
sus, the LMC only contains 12 globular clusters (NGC 1466, 
1754, 1786, 1898, 2210, 2257, 2005, 2019, 1835, 1841, Hodge 
11, and Reticulum, excluding ESO121-SC09, Olsen et al. 
1998; Johnson et al. 1999) with ages comparable to those of 
the oldest Milky Way globulars. This relatively small num- 
ber of objects limits their utility and applicability for the 
study of galaxy-wide properties. 

There is another tool at our disposal for studying the 
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formation and early evolution of the LMC - RR Lyrae vari- 
ables (Smith 1995; Sarajedini 2011). These pulsating stars 
that are located on the horizontal branch are ubiquitous in 
the LMC with tens of thousands having been identified and 
characterized by recent surveys (Alcock et al. 2000; Soszyn- 
ski et al. 2009, see below). They are well known as distance 
indicators, and, in addition, their mere presence necessitates 
the existence of a stellar population with an age greater than 
~10 Gyr. Their periods and amplitudes are useful for mea- 
suring the metal abundance of the star cluster or galaxy in 
which they reside. Furthermore, their minimum light colors 
can be used to measure the line-of-sight reddening. 

There are three principal types of RR Lyrae variables; 
those pulsating in the fundamental mode exhibit sawtooth- 
like light curves and are referred to as ab-type or RRO vari- 
ables. The first overtone pulsators generally show sine-curve 
shaped light curves, have shorter periods and typically lower 
amplitudes than the ab-types, and are referred to as c-type 
or RRl variables. Lastly, RR Lyraes that pulsate in both 
the fundamental and first overtone modes (i.e. double mode 
pulsators) carry the d-type moniker (RRd). 

The most recent survey of RR Lyraes in the LMC is 
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represented by the Optical Gravitational Lensing Experi- 
ment (OGLE, Udalski et al. 1992). We describe the detailed 
properties of this survey in the next section. For now, we 
note that at least four studies have been published utilizing 
the vast number of LMC RR Lyraes from OGLE. Pejcha & 
Stanek (2009) used 9393 ab-type RR Lyraes in the OGLE-III 
catalog to investigate the structure of the LMC stellar halo. 
After adopting a median distance of 50 kpc for their sample 
of RR Lyraes, Pejcha & Stanek (2009) construct a reddening 
map of the LMC and use it to study the line of sight dis- 
tribution of the RR Lyraes. They find that the RR Lyraes 
define a triaxial ellipsoid with axes ratios of 1:2.00:3.50 pro- 
viding further evidence that the LMC does indeed possess a 
halo distribution of RR Lyrae variables (Kinman et al. 1991; 
Minniti et al. 2003; Borissova et al. 2006). 

The series of papers by Haschke et al. (2011), Haschke 
et al. (2012a), and Haschke et al. (2012b) used the OGLE- 
III RR Lyrae data to explore the reddening, metallicity, 
and structure of the LMC, respectively. The first study 
noted above presented reddening maps of the LMC using the 
red clump stars and the RR Lyrae variables and compared 
these maps with several previous studies. The second paper, 
Haschke et al. (2012a), used the Fourier properties of the RR 
Lyrae light curves to study their metallicity distribution in 
the LMC and compared the results with spectroscopy-based 
abundance values from previous studies. The third paper, 
Haschke et al. (2012b), used the line-of-sight distance dis- 
tribution of the RR Lyraes and Cepheids to investigate the 
structural properties of the LMC. We will refer to these stud- 
ies later as we compare our results with theirs. 

In the present work, we seek to study the redden- 
ing, metallicity, and spatial distribution of the OGLE-III 
RR Lyraes using complementary techniques to those em- 
ployed in the previous studies of Pejcha & Stanek (2009) 
and Haschke et al. (2011, 2012a, 2012b). We will compare 
our results with those of previous investigations in order to 
test the veracity of the techniques and the results they yield. 
We also include the relevant properties of the old LMC glob- 
ular clusters in our analysis to gain some insight into their 
relationship with the RR Lyraes. This paper is organized as 
follows. The next section provides a brief description of the 
OGLE-III RR Lyrae database. Section 3 presents our re- 
sults on the metallicity, extinction, and spatial distribution 
of the ab-type RR Lyraes and compares them with previ- 
ous studies. Lastly, Section 4 presents the conclusions of this 
study. 



2 DATA 

The OGLE experiment began its initial observations in 1992 
with OGLE I. Since then, two additional datasets have been 
released (Udalski et al. 1992). The most recent one, from 
OGLE HI, was taken from July 2001 to March 2008 on the 
1.3 meter Warsaw telescope at Las Campanas Observatory 
(Soszynski et al. 2009). The imaging camera is composed 
of eight CCDs each with 2048 by 4096 pixels providing a 
field of view of 35 by 35 arcmin on the sky (Udalski et al. 
2003). Images were taken in the V and I filter passbands. 
The full photometric catalog including the variable objects 
is available online. 

The reduction procedures, photometric calibrations, 



and astrometric transformations are detailed in Udalski et 
al. (2008). The raw light curves of suspected variables were 
subjected to a period search algorithm based on Fourier 
analysis. From this, 24,906 RR Lyrae variables were identi- 
fied and characterized. For each object, the catalog provides 
right ascension, declination, mean magnitudes in I and V, 
period, I-band amplitude, along with the Fourier parameters 
^21, 021, R31; and 031 (Soszynski et al. 2009). 

In the present work, we are mainly interested in the 
17,693 ab-type RR Lyraes. We have excluded those stars 
that were flagged as foreground objects and those with un- 
certain or missing V-band photometry leaving 17,221 ab- 
type RR Lyraes, which constitute the dataset we will ana- 
lyze herein. 



3 RESULTS 
3.1 Metallicities 

Metallicities for our sample are determined using the re- 
lationship from Alcock et al. (2000) shown in equation [TJ 
wherein Pat represents the period of the ab-type variables 
and Vamp is the light curve amplitude in the V-band: 

[Fe/H] = -8.85[log{Pab) + 0.15Kmp] - 2.6. (1) 

To use equation [TJ which is on the Zinn and West (1984) 
abundance scale, we need to convert the I-band amplitudes 
from the OGLE HI catalog to those in the V-band; for this, 
we make use of equation 1 in Dorfi & Feuchtinger (1999). 
They find 

Vamp = 0.075 + 1.497/™p, (2) 

where Vamp and lamp are the light curve amplitudes in the 
V- and I-band, respectively. This empirical relation exhibits 
a R'^ correlation coefficient of 0.904 according to Dorfi & 
Feuchtinger (1999), where R'^ = 1 is a perfect correlation. 

We note that Alcock et al. (2000) quote an error of 0.31 
dex per star in their metallicity-period-amplitude relation 
(equation 1). There is evidence that the actual error may be 
slightly less than this value and closer to 0.25 dex (Jeffery et 
al. 2011). There are at least two additional sources of error 
inherent in equation (1). First, some RR Lyrae stars exhibit 
the Blazhko effect which results in a modulation of their am- 
plitude by as much as 0.2 mag (e.g. Kunder et al. 2010). This 
leads directly to an uncertainty of 0.03 dex in their metal 
abundance. Second, there is the conversion of the I-band am- 
plitude to the V-band value. To estimate the error inherent 
in this, we have examined the original data that were used to 
construct equation (2) from Dorfi & Feuchtinger (1999). We 
have calculated the root-mean-square (rms) deviation of the 
fitted points from the fit (i.e. equation 2) and find a value of 
rms = 0.09 mag in the V band amplitude. Propagating this 
through the Alcock et al. equation yields a contribution of 
0.014 dex to the total metal abundance error. Both of these 
sources of uncertainty are negligible and do not add signif- 
icantly to the uncertainty of 0.31 dex originally quoted by 
Alcock et al. (2000). 

The histogram of the resultant metallicities from the 
Alcock et al. (2000) relation is shown in the left panel of 
Fig. [1] as the solid line. The middle panel of Figure [T] shows 
the metallicity histogram derived using the methodology of 
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Figure 1. The left panel shows the histogram of metallicities 
from the present work (solid line) along with a Gaussian fit to the 
distribution. These are based on the Alcock et al. (2000) equa- 
tion. The middle panel shows the Fourier based metallicities from 
Haschke et al. (2012a) using the Smolec (2005) relation. The right 
panel shows the spectroscopic values from Table 9 of Borrisova et 
al. (2006). The vertical arrows represent the mean metal abun- 
dance of the 12 old globular clusters in the LMC. All of these 
values are on the Zinn & West (1984) scale. 



the Zinn & West (1984) scale and furthermore offset them 
to a common zero point. Since Suntzeff et al. (1992) presents 
metallicities for all of the clusters, we have offset the values 
from the other sources onto the system of Suntzeff et al. 
(1992) and averaged the abundances for clusters in common 
between the various studies. These adjusted values and the 
averages are shown in Table 2. The errors in each average 
value have been calculated using the small sample statistical 
formulae of Keeping (1962). 

The vertical arrow in Fig.[T]shows the mean metal abun- 
dance of the 12 old globular clusters in the LMC taken from 
Table 2. This value equals [Fe/H] = -1.82 ± 0.06 (stan- 
dard error of the mean, sem). We see that this mean value 
is ~0.3 dex more metal-poor than the field RR Lyrae vari- 
ables. However, given the uncertainties in both values, it is 
unclear whether this difference of ~0.3 dex is statistically 
significant. If it is significant, then it could reflect a slight 
age difference between the old globular clusters and the field 
RR Lyrae variables, in the sense that the former would be 
older than the latter. 



Haschke et al. (2012a). Their study utilized the relationship 
derived by Smolec (2005) applied to the Fourier coefficients 
of OGLE RR Lyraes. The Smolec (2005) relation is 

[Fe/H]jg5 = -3.142 -4.902Pa6 -f 0.82,^31, (3) 

but it should be noted that a phase change must be ac- 
counted for in this equation, as the OGLE data uses a cosine 
Fourier series for their light curve fit while Smolec (2005) 
used a sine Fourier series. Equation[3]is on the Jurcsik (1995) 
metallicity scale. This was transformed to the Zinn & West 
1984 scale (ZW). 

[Fe/H]zw = 1.02S[Fe/H]jg5 - 0.242, (4) 

which is taken from the paper by Papadakis et al (2000). 
The histogram of metallicties resulting from equation [4] is 
plotted in the middle panel of Figure [T] as the solid line. 

The right panel of Fig. \T\ illustrates the histogram of 
LMC RR Lyrae metallicities from Table 9 in the spectro- 
scopic study of Borissova et al. (2006). These data were 
transformed from the Harris (1996) scale to the Zinn & 
West scale (ZW) by adding -0.06 dex (Gratton et al. 2004; 
Borissova et al. 2006). 

Fitting a Gaussian function to each metallicity distribu- 
tion in Figure [1] we find the following peak values: -1.70 ± 
0.25, -1.63 ± 0.21, -1.62 ± 0.10, for the present sample, the 
Haschke et al. (2012a) sample, and that of Borrisova et al. 
(2006), respectively. The quoted uncertainties represent the 
1-0" standard deviations in each case. Given the errors, these 
values are in statistical agreement, but it should be noted 
that the Fourier method of Haschke et al. (2012a) yields 
metallicities with a slightly smaller dispersion as compared 
with the Alcock et al. (2000) method that uses periods and 
amplitudes, which we have adopted. 

In order to compare the metallicities of the RR Lyrae 
variables with those of the old globular clusters in the LMC, 
we have assembled abundance values from the literature for 
these 12 clusters as shown in Table 1. In order to combine 
abundance values for the same clusters from different stud- 
ies, we need to ensure that they are on the same abundance 
scale. As such, we have converted all of the metallicities to 



3.2 Metallicity Gradient in the LMC 

It is instructive to examine the radial metallicity variation 
among the stellar populations of the LMC. For each RR 
Lyrae variable in our sample, we calculate its angular dis- 
tance from the center using equation (5), 

cos9 — [sin (90 — 5) * sin (90 — Slmc) * cos (q — olmc)] (5) 
+[cos (90 -5) * cos (90 - Slmc)], 

which is taken from van der Marel & Cioni (2001). In this 
equation, aLMC and Slmc are the right ascension and dec- 
lination of the LMC center from van der Marel & Cioni 
(2001), while a and 5 are the right ascension and declina- 
tion of each OGLE RR Lyrae star. A mean absolute LMC 
distance modulus of 18.55 (see Sec. 3.4) was used to calculate 
the radial location of each RR Lyrae from the center of the 
LMC. The metallicity data were then binned using 0.5 kpc 
bins and plotted in Fig. [5] The binned averages are plotted 
as squares along with a best fit line in Figure (2] This linear 
regression fit to the RR Lyrae data has a slope of -0.027 
± 0.002 dex/kpc, which is steeper than the slope of -0.015 
± 0.003 dex/kpc found in Feast et al. (2010), who used the 
periods of the OGLE RR Lyraes along with equation (1) 
of Sarajedini et al. (2006) in order to calculate metallicities. 
The open circles in Fig.[2]are the old globular clusters in the 
LMC as shown in the last column of Table 2. A least squares 
fit to these globular cluster points yields a slope of -0.022 ± 
0.013 dex/kpc. When the student's T-test is applied to this 
latter value, we find a 93% probability of significance sug- 
gesting that, given the data in Table 2, the radial metallicity 
slope of the LMC globular clusters is not likely to be statis- 
tically significant. We note that the open triangles in Fig. 2 
represent the data from Haschke et al. (2012a) binned in an 
identical manner to our data. The least squares fit to these 
stars yields a slope of -0.019 ± 0.002 dex/kpc, which lies 
between the values from the present study and that of Feast 
et al. (2010). It's not immediately clear why the techniques 
applied by Feast et al. (2010), Haschke et al. (2012a), and 
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Table 1. Original LMC Globular Cluster Metallicity Data 



Cliist6r 


Suntzcff 


Grocliolski 


Mucciarelli 


Colucci 


Na,iiiG 


et al. (1992) 


ct al. (2007) 


et al. (2010) 


et al. (2011) 


NGC 1466 


-1.85 








WGO 1754 


—1.54 








NGC 1786 


-1.87 




-1.75 




NGC 1898 


-1.37 








NGC 2210 


-1.97 




-1.65 




NGC 2257 


-1.80 


-1.59 


-1.95 




NGC 2005 


-1.92 






-1.50 


NGC 2019 


-1.81 


-1.31 




-1.61 


NGC 1835 


-1.79 








NGC 1841 


-2.11 








Hodge 11 


-2.06 


-1.84 






Reticulum 


-1.71 


-1.57 







Table 2. Adjusted LMC Globular Cluster Metallicity Data 



Cluster Suntzeff Grocliolski Mucciarelli Colucci Average 

Name et al. (1992) et al. (2007) et al. (2010) et al. (2011) [Fe/H] 



NGC 1466 


-1.85 








-1.85 


NGC 1754 


-1.54 








-1.54 


NGC 1786 


-1.87 




-1.86 




-1.86 ± 0.01 


NGC 1898 


-1.37 








-1.37 


NGC 2210 


-1.97 




-1.78 




-1.87 ± 0.12 


NGC 2257 


-1.80 


-1.86 


-2.00 




-1.89 ± 0.07 


NGC 2005 


-1.92 






-1.81 


-1.87 ± 0.07 


NGC 2019 


-1.81 


-1.60 




-1.92 


-1.78 ± 0.11 


NGC 1835 


-1.79 








-1.79 


NGC 1841 


-2.11 








-2.11 


Hodge 11 


-2.06 


-2.07 






-2.06 ± 0.01 


Reticulum 


-1.71 


-1.85 






-1.79 ± 0.09 



the present work yield different radial metallicity gradients. 
Further insights on which of these results is closest to reality 
will require further spectroscopic metallicity measurements 
of individual RR Lyraes in the LMC. 

Also plotted in Fig. [2] are the metallicity gradients in 
the disk of the Milky Way and M33 from Friel et al. (2002) 
and Tiede et al. (2004), respectively. The former is based on 
the mean relation for all of the open clusters in the study of 
Friel et al. (2002) while the latter is the mean metallicity of 
RGB stars in an outer disk field of M33. These are meant 
to illustrate a comparison between the disk metallicity gra- 
dients of these three galaxies. In particular, it would appear 
that the LMC relation is significantly shallower than those 
of the Milky Way and M33, which are both very similar. 

One reason for the different metallicity gradients in the 
disks of these galaxies could be related to the ages of the 
stellar populations used to measure these gradients. Friel et 
al. (2002) found that older open clusters in the Milky Way's 
disk exhibit a significantly steeper radial metallicity gradi- 
ent as compared with younger clusters. This trend seems 
to be reversed in the present case where the oldest stars in 
the LMC have a shallower metallicity gradient as compared 
with generally younger populations in the Milky Way (open 
clusters) and M33 (RGB stars). 

Another factor that could influence the slope of the ra- 
dial metallicity gradient in a galaxy is the presence or ab- 
sence of a stellar bar. As Zaritsky et al. (1994, see also Gro- 



cholski et al. 2007, Sharma et al. 2010) point out, barred spi- 
ral galaxies are more likely to possess shallower metallicity 
gradients as compared with non-barred spirals. As such, this 
would explain the appearance of Fig. 2, which suggests that 
the strong bar structure of the LMC has acted to diminish 
the extent of its disk metallicity gradient as compared with 
M33, which has no measurable stellar bar, and the Milky 
Way, which possesses a relatively small stellar bar. 



3.3 Extinction Mapping 

In order to determine the reddening for each RR Lyrae in 
our sample, we make use of the results of Guldenschuh et al. 
(2005) who found that the minimum light color of ab-type 
RR Lyrae variables is equal to {V — I)o ~ 0.58 ± 0.02 regard- 
less of their intrinsic properties. As a result, we can calculate 
the reddening of an individual RR Lyrae via equation [6] as 

E{V -I) = {V- /),„„, - 0.58. (6) 

To estimate the minimum light colors of the ab-type RR 
Lyraes in the OGLE III database, we proceed as follows. 
If the light curves were symmetric sine-like curves, we could 
calculate the minimum light color via the following equation: 

(V - /)_ = [< V > - [< / > +^]. (7) 
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Table 3. LMC Globular Cluster Reddenings 



Cluster E(B-V) E(V-I) Reference 



NGC 1466 





09 





12 


Johnson et al. (1999) 


NGC 1754 





09 





12 


Olsen et al. (1998) 


NGC 1786 





07 





09 


Brocato et al. (1996) 


NGC 1898 





07 





09 


Olsen et al. (1998) 


NGC 2210 





06 





08 


Brocato et al. (1996) 


NGC 2257 





04 





05 


Johnson et al. (1999) 


NGC 2005 





10 





13 


Olsen et al. (1998) 


NGC 2019 





06 





08 


Olsen et al. (1998) 


NGC 1835 





08 





10 


Olsen et al. (1998) 


NGC 1841 





18 





23 


Brocato et al. (1996) 


Hodge 11 





08 





10 


Johnson et al. (1999) 


Reticulum 





05 





07 


Mackey & Gilmore (2004) 
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Figure 2. The radial metallicity gradient of the LMC is shown 
wherein the open squares represent the binned radial distribution 
of the OGLE RR Lyraes in the present study while the open cir- 
cles are the old LMC globular clusters from Table 2. The solid 
line is the least-squares fit to the RR Lyrae relation while the dot- 
dashed line is the fit to the globular clusters. The triangles repre- 
sent the RR Lyrae abundance data from Hashcke et al. (2012a) 
binned in an identical manner to our metallicity data with the 
dotted line representing a least squares fit to these points. The 
metallicity gradients of M33 and the Milky Way are indicated 
by the dotted and dashed lines, respectively. The former is taken 
from the work of Tiede et al. (2004) while the latter represents 
the mean relation for the open clusters in the work of Friel et al. 
(2002). 



However, because the shapes of ab-type RR Lyrae light 
curves are asymmetric, we need to determine the correc- 
tion to be applied to the minimum light colors derived using 
equation (7). In order to calculate this correction, we have 
used the ab-type RR Lyrae light curve templates of Layden 
& Sarajedini (2000) to calculate the quantity 



The mean value of A for the 6 ab-type templates is (A) — 
-0.061 ± 0.017. As a result, we can correct the minimum 
light colors derived from equation (7) by adding -0.061 mag 
to each of our values. 

This corrected apparent minimum light color is then 
coupled with the results of Guldenschuh et al. (2005) shown 
in equation (6) above to calculate the line-of-sight redden- 
ing. Figure ^illustrates the resultant reddening distribution 
for the LMC RR Lyraes in our sample. The number of RR 
Lyraes exhibiting a negative reddening in Fig. |3] represents 
only 5% of the total, which is a testament to the robustness 
of our method (Sarajedini et al. 2006). The error in each 
reddening is composed of a— 0.02 mag from the uncertainty 
in the value of {V — I)o,min and (t=0.017 mag from the cor- 
rection applied to the minimum light colors. Taken together, 
these combine to produce an uncertainty of (t=0.026 mag in 
E(V-I). 

In order to test our method for reddening estimation, 
we have extracted the V and I band time-series photomet- 
ric data for 200 RR Lyraes from the OGLE III survey. We 
have performed light curve fitting on these data (Layden & 
Sarajedini 2000) using the known periods in order to derive 
the minimum light colors of these stars. These were then 
used to calculate reddenings using equation 6 above. When 
we compare these reddenings to the results of applying our 
method (equations 7 and 8), we find a mean difference of 
0.0291 ±0.005 (sem) mag in E(V-I). This is consistent with 
the error we have quoted for our reddening determinations 
of ±0.026 mag thus providing support for the technique we 
have used to calculate the line of sight reddening to the LMC 
RR Lyraes. 

A Gaussian profile was fit to the distribution of redden- 
ing values for the RR Lyrae. The peak of this distribution 
is at E(V-I) = 0.12 with a = 0.05, which corresponds to 
E(B-V) = 0.09. This value is in good accord with the work 
of Haschke et al. (2011), who found a mean LMC reddening 
of E(V - I ) = 0.09 ± 0.07 mag from their analysis of red 
clump stars in the LMC and E(V - I ) = 0.11 ± 0.06 mag 
based on the OGLE RR Lyrae stars. In fact, we can perform 
a star-by-star comparison of our reddenings with those from 
Haschke et al. (2011). Doing this, we find a mean differ- 
ence of -0.021 ± 0.041, which suggests that the two datasets 
are statistically indistinguishable. This offset is also consis- 
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E(V-I) 

Figure 3. Histogram of reddening values for LMC ab-type RR 
Lyraes in our data set. The solid line is the Gaussian fit to the 
central portion of the distribution. 

tent with the quoted error of ±0.026 in our values of E(V- 
I). It should be noted that Haschke et al. (2011) utilized a 
completely different method than we have to estimate the 
RR Lyrae reddenings. In particular, they rely on RR Lyrae 
metallicities derived from the Fourier coefficients, which are 
then used in conjunction with theoretical light curve models 
to infer the intrinsic color of each RR Lyrae star in their 
sample. 

It is also interesting to note that the mean reddening of 
the 12 old globular clusters in the LMC is E(V-I) = 0.11 ± 
0.06 (see Table 3), which is consistent with that of the RR 
Lyraes. All of these reddening values are, in turn, consistent 
with the range of values published in Table 2 of Imara & 
Blitz (2007), which extends from E(V-I) « 0.03 to E(V-I) « 
0.26 with a mean value of E(V-I) « 0.16. 

With reddenings for over 17,000 ab-type RR Lyraes in 
hand, we can also investigate the extinction map of the 
LMC. This is shown in Fig. |4l wherein the squares rep- 
resent the locations of the old globular clusters in the LMC 
from Table 1 of Soszynski et al. (2009) . There are a number 
of features of note in this diagram. First, the area around 30 
Doradus and the southern CO ridge are prominent sources 
of extinction in the case of the RR Lyraes. There is an ad- 
ditional extinction ridge extending north-south along the 
western edge of the LMC surveyed area at RA«4''50™. It 
is also important to note that there is an extinction "hole" 
in the central region of the LMC. Generally speaking, the 
oldest globular clusters in the LMC (blue squares in Fig. [3)) 
tend to lie in areas of reduced extinction. This phenomenon 
brings forth the question of whether undiscovered old clus- 
ters may exist in areas where high extinction has prevented 
their identification. If so, it is vitally important for these old 
clusters to be identified and characterized because the over- 
all dearth of old clusters in the LMC has made it challenging 
to study its early formation epochs using these clusters. 

With a reddening map for the LMC in hand, it is of 
interest to compare it with maps of the CO emission from 
the LMC. The latest version of such a map is from the work 
of Wong et al. (2011) representing results of the Magellanic 
MOPRA Assessment (MAGMA) survey. This survey builds 
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X (degrees from center) 

Figure 4. The contour map of extinction across the LMC based 
on the minimum light colors of the ab-type RR Lyraes from the 
OGLE III survey. The small points indicate the locations of the 
individual RR Lyrae variables while the blue squares are the old 
globular clusters in the LMC listed in Table 3 of Soszynski et al. 
(2009). 




X (degrees from center) 

Figure 5. The contour map of LMC CO molecular cloud inten- 
sites (K km/s) from the MAGMA survey of Wong et al. (2011). 
These are plotted on the same scale as the RR Lyrae reddening 
map. 

upon the previous work by Fukui et al. (2008) who used the 
4.0m NANTEN radio telescope. Fig. [5]shows the intensity in 
units of K km/s of the CO emission from the clouds studied 
by Wong et al. (2011) from their Table 6. This figure is 
plotted to the same scale as Fig. |4l which represents the 
RR Lyrae reddening map. 

A casual comparison of Fig. [5] and Fig. |4] suggests a 
high degree of correlation between the extinction as revealed 
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3000 



X (degrees from center) 

Figure 6. The RR Lyrae reddening map of tlie LMC is stiown in 
color. Overplotted in greyscale is the contour map of CO molec- 
ular cloud intensities (K km/s) from the work of Wong ct al. 
(2011). 



by the RR Lyrae variables and the intensity of CO present 
in the line of sight. Further insight into this correlation can 
be gained by overplotting the two distributions as shown 
in Fig. [6l Once again, this figure emphasizes that the RR 
Lyrae reddenings and the CO intensity distributions are well 
correlated in their on-sky distribution. This suggests that 
the colors of the RR Lyraes are influenced by dust in the line 
of sight further indicating that some of the RR Lyraes are 
located behind and within the CO dust clouds. The one area 
that seems to be an exception to this assertion is the central 
region of the LMC, where we see emission from CO clouds 
but where the RR Lyrae reddenings are all relatively small 
(see Fig. It is possible that, in this area, the CO clouds 
are predominantly on the far side of the LMC and/or the 
RR Lyraes are on the near side so that the dust extinction 
does not influence the colors of the RR Lyrae stars. 



3.4 Distances 

With metallicities, reddenings, and mean magnitudes in 
hand, it is a simple matter to calculate the distance to each 
of the over 17,000 ab-type RR Lyraes in the OGLE sam- 
ple. First, we estimate the absolute magnitude of each RR 
Lyrae using = (0.23±0.04)([Fe/i/] + 1.6) + (0.56±0.12) 
from Chaboyer (1999), and calculate the extinction using 
Av = 2.38£(V - /) from Schlegel et al. (1998). 

The resulting distribution of RR Lyrae distance mod- 
uli is shown as the open circles in Fig. [T] A Gaussian 
proflle fit to this distance distribution yields a peak value 
of (m — M)o = 18.55 with a 1-a standard deviation of 
0.10 mag. This translates to a distance of 51.3 ± 2.4 kpc. 
This value is well within the range of values for the dis- 
tance to the LMC presented in the compilation of Clemen- 
tini et al. (2003). Alternatively, we can perform this cal- 
culation using the Benedict et al. (2011) calibration be- 



2000 



1000 



■ 




li 












\ B 








7 * 


t 





18.4 



18.8 



(m-M)o 



Figure 7. The open circles represent the histogram of distance 
moduli for the RR Lyraes in our sample using the Chaboyer et 
al. (1999) relation between RR Lyrae absolute magnitude and 
metallicity. The solid line is the Gaussian distribution fitted to 
the central portion of this histogram. The vertical dashed line 
represents the peak of the fitted Gaussian distribution while the 
dotted lines show the it 1 kpc range of the putative LMC disk pop- 
ulation (see text). In comparison, the filled circles represent the 
histogram of distance moduli yielded by the Benedict et al. (2011) 
relation between RR Lyrae absolute magnitude and metallicity. 



tween RR Lyrae absolute magnitude and metal abundance 
- Mv = (0.214 ± 0M7){[Fe/H] + 1.5) + (0.45 ± 0.05). This 
yields the binned histogram shown by the filled circles in 
Fig. [T] which has a peak value of (m — M)o = 18.68 with 
a 1-a standard deviation of 0.10 mag, as yielded by a Gaus- 
sian fit. Given the errors, this distance modulus is consistent 
with the one derived from the Chaboyer (1999) RR Lyrae 
relation, but it is on the high end of the distribution of LMC 
distances from Clementini et al. (2003). 

The question arises as to whether the line-of-sight dis- 
persion in distance is statistically significant or simply the 
result of measurement errors. To address this question, we 
note that Alcock et al. (2000) estimate metallicity errors of 
0.31 dex from equation (1) above. This translates to an er- 
ror of (7=0.071 mag in the absolute magnitude of each RR 
Lyrae. However, it is important to keep in mind that the 
1-a dispersion of the RR Lyrae metallicity distribution de- 
rived from the Alcock et al. (2000) relation (Fig. 1) is 0.25 
dex, which is smaller than the error of an individual mea- 
surement of 0.31 dex. This suggests that 0.31 dex may be 
an overestimate of the error in equation (1). In any case, 
as discussed above, the reddening error is equal to ±0.026 
mag in E(V-I) leading to an uncertainty of (t=0.062 mag 
in the determination of the interstellar extinction. Ignoring 
for the moment the uncertainties in the measured V-band 
magnitudes, we find a distance modulus error of ct=0.094 
mag. This error is consistent with the dispersion measured 
in the distance histogram shown in Fig. [7] suggesting that 
this line of sight dispersion in distance is dominated by the 
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measurement errors in metallicity and reddening. We return 
to this point in the next section. 

3.5 Distance Along Line of Maximum Gradient 

To examine the structure of the LMC as revealed by the 
ab-type RR Lyraes, we can use the "kinematic circular disk 
method" as detailed in van der Marel & Cioni (2001) and 
revisited by Grocholski et al. (2007). In this method, one 
assumes a model of a circular disk with stars in circular and 
non-spherical orbits. Then, when plotting the line-of-sight 
distance versus the distance along the line of the maximum 
velocity gradient, objects following circular orbits should fall 
along a distinct line. The slope of this line should indicate 
the inclination of the disk. Thus, objects which deviate from 
this slope likely have non-circular orbits. 

In Figure (S] we show our derived distances for the RR 
Lyrae stars as a function of the distance along the line of 
maximum velocity gradient. Also plotted are the clusters 
from Grocholski et al. (2007) and Walker (1992). Figure [H] 
is modeled after Fig. 5 of Grocholski et al. (2007), except 
that here we plot the line-of-sight distances relative to the 
mean distance of the LMC. While the vast majority of the 
clusters in Figure [8] are found to reside within ~1 kpc of 
the LMC disk (as indicated by the dashed lines), the RR 
Lyrae variables do not appear to follow this trend. The dis- 
tribution of RR Lyraes shows no indication of being oriented 
along the slope of the line (i.e. the inclination of the LMC 
disk). In addition, a significant number of the RR Lyraes 
(f»45%) are found in the region beyond the 1 kpc thick disk 
(see Fig. Q assumed for the LMC (Grocholski et al. 2007). 
One interpretation of this result is that the RR Lyraes in the 
OGLE sample populate a kinematically hot halo around the 
LMC, consistent with the findings of Minniti et al. (2003) 
and Pejcha & Stanek (2009, for a contrasting view, see Sub- 
ramaniam & Subramanian 2009). However, as discussed in 
the previous section, the line-of-sight distance distribution 
of the RR Lyrae seems to be dominated by the measurement 
errors. As a result, the value of ~45% is likely to be an up- 
per limit for the fraction of the OGLE ab-type RR Lyraes 
that are members of the LMC halo. 



4 CONCLUSIONS 

In this paper, we present metallicities, reddenings, and dis- 
tances of RR Lyrae stars in the LMC based on the OGLE 
III database. These stars are used to examine the metal- 
licity gradient of the LMC, its extinction profile, and the 
line-of-sight distribution of the RR Lyrae as a proxy for the 
LMC's halo population. From this investigation, we draw 
the following conclusions: 

1. There appears to be a radial metallicity gradient of 
the RR Lyraes in the LMC. The significance of the gradient 
depends on how the properties of the RR Lyrae light curves 
(e.g. periods, amplitudes, Fourier coefficients) are converted 
to metallicity. Given this apparent dependence, it is impor- 
tant to measure spectroscopic metal abundances for a large 
sample of LMC RR Lyraes with a substantial radial extent. 

2. The extinction of the RR Lyraes in the LMC is cor- 
related with the location and emission of the CO clouds 
detailed by previous studies. However, the CO clouds do 
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Figure 8. Line of sight distance versus distance along the line 
of maximum velocity gradient for a circular disk model (both 
in kpc). RR Lyraes from the present work are plotted as dots; 
old globular clusters are shown as diamonds (Grocholski et al. 
2007) and triangles (Walker 1992). The solid line represents the 
expected slope of the LMC center. The dashed lines represents 
the 1 kpc thick disk of the LMC. 

not appear to be correlated with the density of RR Lyraes 
across the face of the LMC. In addition, the old LMC glob- 
ulars tend to be located in regions of low extinction. This 
underscores the need to survey the regions of higher extinc- 
tion searching specifically for more old globular clusters. 

3. Our analysis of the OGLE III ab-type RR Lyraes 
is inconclusive with regard to whether these stars belong 
to the halo or disk of the LMC. The RR Lyraes appear to 
have a line-of-sight distribution that is not consistent with 
the expected behavior of the LMC disk. However, the un- 
certainties inherent in our analysis and the distribution of 
the OGLE RR Lyraes in the LMC render this a tentative 
conclusion. In any case, previous work, such as Minniti et al. 
(2003), does suggest that there are RR Lyraes in the LMC 
that belong to a kinematically hot halo. 
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